Variability in GRB Afterglows and GRB 021004 



Ehud Nakar a>1 , Tsvi Piran a ' , Jonathan Granot ' 

A Racah Institute for Physics, The Hebrew University, Jerusalem, 91904, ISRAEL 
b Institute for Advanced Studies, Princeton, NJ 08540, USA 



Abstract 

We present general analytic expressions for GRB afterglow light curves arising from 
a variable external density profile and/or a variable energy in the blast wave. The 
former could arise from a clumpy ISM or a variable stellar wind; The latter could 
arise from refreshed shocks or from an angular dependent jet structure (patchy 
shell). Both scenarios would lead to a variable light curve. Our formalism enables 
us to invert the observed light curve and obtain possible density or energy profiles. 
The optical afterglow of GRB 021004 was detected 537 seconds AB (after the burst) 
(Fox et al. 2002). Extensive follow up observations revealed a significant temporal 
variability. We apply our formalism to the R-band light curve of GRB 021004 and 
we find that several models provide a good fit to the data. We consider the patchy 
shell model with p = 2.2 as the most likely explanation. According to this model 
our line of sight was towards a "cold spot" that has lead to a relativity low 7-ray 
flux and an initially weak afterglow (while the X-ray afterglow flux after a day was 
above average). Observations above the cooling frequency, v c , could provide the best 
way to distinguish between our different models. 
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1 Introduction 



The behavior of gamma-ray burst (GRB) afterglows is well known for a spher- 
ical shell propagating into a constant density inter-stellar medium (ISM) or 
into a circum-burst wind with a regularly decreasing density. Sari, Piran & 
Narayan (1998, hereafter SPN98) have presented a simple analytic model for 
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the ISM case, assuming synchrotron emission from an adiabatic relativistic 
blast wave. Chevalier & Li (1999) generalized this model for a circum-burst 
wind density profile. In both cases the flux shows a spectral and temporal seg- 
mented power law behavior, F v oc t a i/P. The indices a and (3 change when a 
spectral break frequency (the cooling frequency, u c , the synchrotron frequency, 
v m , or the self absorption frequency, u sa ) passes through the observed band. 
The values of the spectral and temporal indices depends on the cooling regime 
(fast or slow) and on the ordering of v relative to v sa , v c and v m . Most GRB 
afterglows display a smooth power law decay. 

In several cases the observed afterglow light curves have shown deviations 
from a smooth power law. The most prominent case is the recent GRB 021004 
whose optical counterpart was observed at a very early time, 537 sec after 
the trigger (Fox et al. 2002). Following observations at short intervals showed 
fluctuations around a power law decay. We develop here the general theory 
for GRB afterglows when the relativistic blast wave encounters a variable 
external density or its energy (per unit solid angle) varies with time. Such 
variations in energy could arise due to refreshed shocks, when initially slower 
moving matter encounter the blast waves after it has slowed down (Kumar & 
Piran, 2000a), or due to angular variability within the relativistic jet (Kumar 
&: Piran, 2000b). Both variation in the density or in the energy can reproduce 
a variable light curve and in particular the observed R-band light curve of 
GRB 021004. However, as we argue latter, there are some weak indications 
that a variable energy model that arises from a patchy shell structure (random 
angular fluctuations in the jet) seems to give the best fit to all the available 
data. If correct this interpretation implies that the electron power law index 
is p ~ 2.2, a suggestion that might be confirmed with a more detailed multi- 
wavelength spectrum. 



2 Theory 



We generalize the results of SPN98 to a time dependent energy and a spa- 
tially varying external density. We first outline the general model and then 
investigate two specific cases. Following SPN98 we assume that the dominant 
radiation process is synchrotron emission. In our model, the mass in the blast- 
wave at radius R is taken to be the integrated external mass up to this radius, 
and we assume that all this mass is radiating. The internal energy density 
of the emitting matter at radius R is taken from the shock jump conditions, 
which depend only on j(R) and n(R). These approximations are valid as long 
as the external density and the energy in the blast-wave do not vary too 
rapidly. For example a large density jump can produce a reverse shock while 
a sharp density drop may initiate a rarefaction wave. The accuracy of this 
model decreases as the variations in the density and the energy become more 
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rapid. 



A few hundred seconds after the GRB the relativistic ejecta decelerates, driv- 
ing a strong relativistic shock into the ambient medium. As radiative losses 
become negligible the flow settles into the adiabatic self-similar Blandford- 
McKee (1976) blast wave solution. Energy conversion takes place within the 
shock that propagates into the external medium. The energy equation reads: 

E(t) = A^[R(t)]M[R(t)}c 2 , (1) 



where E is the isotropic equivalent energy and A is a constant of order unity 
whose exact value depends on the density profile behind the shock (e.g. for 
an external density n(r) oc r~ k , A = 4(3 — fc)/(17 — 4k); Blandford & McKee 
1976). In the following we use A = 1. In equation (1), M(R) is the mass of 
the blast-wave, i.e. the integrated external mass up to a radius R, 

R 

M(R) = An J n(r)r 2 dr . (2) 
o 



The observed time, t, is related to R and 7 through two effects. First, the 
observed time of a photon emitted on the line of sight at a radius R is t\ os = 
(l/4c) / H 7~ 2 dr. Second, photons emitted at different angles at the same radius 
R are observed during an interval of ~ R/2c^ 2 . Following SPN98 we estimate 
the observed time interval during which most of the emission emitted at radius 
R is received as t ang ~ R/Acy 2 . Therefore: 



(3) 




For a constant density ISM t ang = 4ti os and t = 5R/16cy 2 w RjArfc. Of 
course, this treatment of the angular effects is only approximate 4 . In most 
cases angular spreading will smooth out any variability on time scales shorter 
than R/Acy 2 . 

We further assume that the electron energy distribution is a power law with an 
index p, and that the magnetic field and the electrons hold fractions 6b and e e , 
respectively, of the internal energy. Now, taking 1^(7, n, €b, e e ), ^(7, n, e#, t) 
and F v ^ max {M,~f,n,€B) from SPN98 and the equations above we obtain: 

v m = 5 • l^E^M^nl^eT.A-, Hz > ( 4 ) 

4 See the Appendix for an extended discussion of the angular smoothing effect and 
Nakar & Piran (2003) for a solution that takes a full account of this effect 
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v c = 3 • 10 u E^Mi 9 n^ /2 tfes% 2 Hz , (5) 
i^max = 7 E 52 nl /2 e)i 2 2 D 2 £ mJy , (6) 



where Q x denotes the value of the quantity Q in units of 10 x (c.g.s), td is the 
observed time in days, D is the distance to the GRB, and for simplicity we do 
not include cosmological effects throughout the paper. The above equations 
readily provide expressions for the flux density at different frequencies 



u (l-p)/2 E P M l-P n (l+p)/4 e (l+PV± e (p-l) Um <U<U c 

F v (x I , (7) 



We concentrate on the above two power law segments, since they are usually 
expected to be the most relevant for the optical light curve. Similar expressions 
for other power law segments of the spectrum may be derived similarly. 

These are the generic expressions for a varying energy and a varying external 
density profile. In addition to the explicit dependence on t in Eq. 7 there is 
an implicit dependence through E(t), M[R(t)] and n[R(t)}. For an ISM or 
wind, M[R(t)] and n[R(t)} have simple analytic forms and Eq. 7 reduces to 
the expressions of SPN98 and Chevalier k Li (1999). 

For p ~ 2, F v>Vc is only weakly dependent on M and n, while the dependence 
on E is roughly linear (note that F u depends on E also implicitly through R 
that appears in M(R)). This feature enables us to distinguish between energy 
dominated fluctuations and density dominated fluctuations in the afterglow 
light curve, when there are measurements both above and below the cooling 
frequency, v c . 

In reality, it is unlikely that both variations (in E and in n) will be important 
in a given burst (since this would require a coincidence). Therefore, we shall 
consider below, in some detail, the cases where one of these quantities is con- 
stant while the other one varies. Moreover, the information in a single band 
light curve (or more accurately, from a single power law segment of the spec- 
trum) is insufficient to determine both profiles. For any given set of density 
and energy profiles the light curve can be easily calculated. However, these 
profiles are not at hand. The observable is the light curve and these profiles 
are unknown variables. It is necessary to make some assumption for one of 
the profiles in order to deduce the other (for example, to assume a constant 
energy or a constant density). 
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2.1 A Variable External Density 



Consider, first, the case where the dominant variations are in the density 
profile while the energy is constant. Eqs. 2, 3 and 7 reduce to: 



F v oc < 



M l-p n (l+p)/4 m< < 

(8) 

M 2-p n (p-2)/4 r < 



1 



J Mdr , 
o / 



t = ^\MR+ I Mdr] , (9) 
M = Anm p J r 2 n(r)dr . (10) 



For a given -F„(t) we solve Eqs. 8-10 for R(t), n[R(t)] and M[R(t)] with pas a 
free parameter. The integral dependence of M[n(r)] in Eqs. 9 and 10 makes it 
difficult to invert these equations analytically for an arbitrary density profile 
(an exact numerical solution is always possible). However, an approximate 
analytic solution can be obtained if the density profile varies slowly (note that 
as discussed earlier, when the density varies rapidly our whole approach is less 
accurate). 

As M grows monotonically with R, t ang is always larger than t\ os and we can 
approximate t ~ t ang . Taking the time derivative of Eq. 8 for v m < v < u c and 
using Eq. 10 we obtain: 

d\nF u (l+p)d\nn 

—77— = m ^ -. -77— for V m < V < V c , (11) 

dint 4 dint 



where 5 m = (1 —p)/[l + n/(3n)] and n(R) = M/(4/3)7im p R 3 is the average 
initial density inside a sphere of radius R. If 5 m varies slowly with time we 
derive: 

F v = F {t/t ) 5m (n/n ) (1+p)/4 for u m < v < v c , (12) 



where F and no are the flux and density at some given time to. As long as 
n > n, 5 m depends weakly on n/n and its value varies between 1 — p and 
0.75(1 — p). When n < n, 5 m ^ and F v depends on t only implicitly via n. 

A similar derivation for u c < v results in: 

F v = F (t/t ) 5c (n/n ) (p ~ 2)/4 for v c <v, (13) 
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where 5 C = (2— p) /[l+(l/3)n/n] — 1. The explicit dependence onn is negligible, 
n O~ 2 )/ 4 ; an( j the variations in n/n yield 1— p < 5 C < — 1. The variations in 5 C 
could be measured if p is large and n 3> n. However, in this limit the density 
changes very rapidly, so that our formalism may not hold. Both Eqs. 12 and 
13 contain the wind solution (with 3n = n and n oc R~ 2 oc t~ r ) and the ISM 
solution (with n = n ). 



2.2 A Varying Energy 



Consider now the afterglow when the energy in the emitting region varies 
with time but the density profile is regular. In the ISM case Eqs. 7 and 3 are 
reduced to: 



F u oc ISM , (14) 

E P-i R 3(2- P ) t -i Uc<v 

' = \e + / E dr ) ' (15 » 



In the wind case ( n = A W R ) these equations become: 

FvX \ i i«2 ^ i (wind) ' (16) 

t = nA w m p c ^ + / ^ d ^j ■ ( 17 ) 



Again, these equations can be solved numerically for a given F u (t). Note that 
in this case the condition t ang > t\ os does not always hold. A sharp increase in 
E would decrease t ang without affecting t\ os . However, if the energy profile is 
not too steep, the condition t ang > t\ os does hold, and we can approximate t 
by t ang . In this case, Eqs. 14 and 16 reduce to the well known ISM and wind 
equations for a constant energy, where E is replaced by E(t). 

Two different phenomena could cause energy variations in the emitting region: 
refreshed shocks and initial energy inhomogeneities in the jet. Refreshed shocks 
(Kumar & Piran 2000a) are produced by massive and slow shells, ejected late 
in the GRB, that take over the blast-wave at late times, when the blast- 
wave has decelerated. These shells bring new energy into the blast-wave. The 
collision produces a refreshed forward shock propagating into the blast-wave 
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and a reverse shock propagating into the slower shell. After these shocks cross 
the shells the blast-wave relaxes back to a Blandford & McKee (1976) self- 
similar solution with a larger total energy (Since the mass of the blast wave 
is dominated by the swept circum-burst material, we neglect the mass of the 
inner shell). At this stage the observed flux is similar to the one emitted 
by a constant energy blast wave with the new and larger energy. Refreshed 
shocks can only increase the energy. Therefore a refreshed shocks energy profile 
should grow monotonically with time, most likely in a step wise profile (each 
step corresponds to the arrival of a new shell). 

Initial energy inhomogeneities (the patchy shell model of Kumar & Piran 
2000b) in the jet could be either regular or irregular ones. During the jet evo- 
lution regions within the relativistic flow with an angular separation larger 
than 7 _1 are casually disconnected. Therefore, the inhomogeneities could be 
smoothed only up to an angular scale of 7 _1 . As 7 decrease the causal con- 
nected regions grow and the initial inhomogeneities can be smoothed on an- 
gular scale of 7 _1 . Recent numerical hydrodynamical studies (Kumar & Gra- 
not 2002) show that at early times the initial fluctuations remain almost un- 
changed, and are smoothed only at rather late times. Additionally, due to 
relativistic beaming, an observer can see only regions within an angle of 7 -1 
around the line of sight. However, regardless of the degree of hydrodynami- 
cal smoothing of the initial fluctuations, when combined with the relativistic 
beaming, the two effect cause F v (t) to reflect the initial physical conditions 
within a solid angle of ~ 7~ 2 (t). As a consequence, the average energy in 
the observed area varies with 7 and therefore with t. This behavior can be 
approximated by the solution presented above, where E(t) is the averaged 
initial isotropic equivalent energy within a solid angle of 7~ 2 , E(t). 

In the patchy shell scenario, fluctuations would appear in the energy pro- 
file when 7 _1 increases to the typical angular size, 9r, of the initial inhomo- 
geneities. When 7" 1 ~ 6& the nearest neighboring fluctuations begin to be 
observed, and the amplitude of the fluctuations in E (and correspondingly in 
F v ) are largest, of the order of the amplitude of the individual fluctuations, 
Aq. As 7 decreases below O^ 1 , the observed number of fluctuations becomes 
large, Nq ~ (7^h) -2 , and the amplitude of the fluctuations in E decreases to 
~ AqNq ~ A fl 7# fl oc 7. For v > min(z/ m , v c ), F v has a close to linear depen- 
dence on E(t) ~ E(t), so that the amplitude of the fluctuations in F u should 
be similar to those in E, with only minor differences between the different 
power law segments of the spectrum. 

A single bump in the light curve can be seen for an axially symmetric struc- 
tured jet, by an observer at an angle o b s from the jet symmetry axis, at the 
time when 7 -1 ~ # b s . At this time the brighter portion of the jet, near its 
symmetry axis where the energy per unit solid angle is largest, becomes visible 
to an observer at angle # bs- Additional bumps are more difficult to produce. 
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3 The Light Curve of GRB 021004 



GRB 021004 is a faint long burst detected by Hete-2 Fregate instrument. 
The burst redshift is z=2.232 (Chornock & Filippenko, 2002) and its isotropic 
equivalent energy is 6 • 10 52 ergs (Lamb et al. 2002 and Malesani et al. 2002). 
An optical counterpart was first observed 537 sec AB (after the burst) (Fox 
et al. 2002) at an R magnitude of 15.5. After a short power law decay, at 
t ~ 2000 sec, a clear bump (about 1.5 mag above the power law decay) is 
observed. >From this time on, frequent observations showed a fluctuating 
light curve (possibly above and below a power law decay). The inset of Fig. 
1 shows the R-band light curve up to 5 days after the trigger. Observations 
after 6 days show a steepening of the light curve which may be interpreted 
as a jet break (Malesani et al. 2002). A break at this time implies a total 
energy (after beaming corrections) of 5 3 • 10 50 ergs. Chandra observed the X- 
ray counterpart of GRB 021004 at 20.5 hr AB for a duration of 87 ks (Sako k 
Harrison 2002). The corresponding mean 2-10 keV X-ray flux in the observer 
frame is 4.3 • 10~ 13 erg cm -2 sec -1 . The X-ray observations showed a power 
law decay index of —1 ± 0.2 and a photon index of 2.1 ± 0.1 which imply an 
electron index p = 2.2 ± 0.2. 

We use the two models described above to find a varying density profile or a 
varying energy profile that reproduce the light curve of GRB 021004. We fit 
the R-band light curve that has the most detailed data. Unfortunately, the 
data in the other bands is not detailed enough and the effect of reddening 
is unknown so a multi wavelength fit is impossible at this stage. We assume 
that the R band is above the synchrotron frequency, z/ m , and below the cooling 
frequency, v c . This assumption is marginal at the time of the first bump (Both 
the transition from fast to slow cooling and the passage of v m through the 
optical bands occur approximately at this time). However, this assumption is 
certainly valid during the later fluctuations of the light curve 6 . It has been 
suggested that v c passes through the optical at t ~ 1 — 3 days (Matheson et al. 
2002). In this case, we expect the fluctuations in the light curve to decrease 
dramatically at t > 3 days, if they are due to fluctuations in the external 
density. We discuss only variability above a constant ISM density profile. As 
we show latter, a reasonable fit with a background wind profile requires an 
electrons' index p < 2 (for either variable density or variable energy), which 
we consider to be a not very physical value. 



5 This value is obtained using a redshift of 2.323 and isotropic equivalent energy of 

6 • 10 52 ergs. The rest of the parameters are similar to those of Frail et. al. 2001. 

6 It is possible that the origin of the first bump is different from the later fluctuations 
(e.g. a passage of v m through the R band combined with the emission from the reverse 
shock, Kobayashi & Zhang 2002), but following Occam's razor we are looking for a 
single explanation to the whole light curve. 
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3.1 A Variable Density Profile 



Lazzati et al. (2002) suggest that the fluctuations seen in the R-band light 
curve arise from variations in the external density profile. They calculate nu- 
merically the resulting light curve for a given density profile, assuming p — 2, 
and show that it agrees with the observations. We invert the observed R-band 
light curve, both analytically and numerically, and derive several possible den- 
sity profiles for different values of p. 

We begin the fit at the first observation, t = 537 sec after the trigger, and 
define no and Rq as the values at this time. For simplicity, we assume a constant 
density up to Rq [so that n(R < Rq) = n ]. With this assumption the ratios 
R/Ro and n/rio do not depends on the values of Rq and no- Figure 1 depicts 
the density profile for a few values of the electron power law index, p. The thick 
lines show the exact numerical solution of Eqs. 8-10, while the thin lines show 
the analytic solution of Eq. 12 (In this solution the value of 5 m is recalculated 
every time step). In order to reproduce the light curve with p > 2.4, the density 
profile must increase with R almost monotonously. Such a density profile does 
not look feasible. For p = 2.2 the density increases by an order of magnitude 
at R ~ 1.5.Ro ~ 3x 10 17 cm and remains roughly constant at larger radii. This 
is consistent with the termination shock of a stellar wind that interacts with 
the ambient medium (Wijers 2001), provided that the latter has a very high 
density of ~ 10 4-5 cm -3 in order for the radius of the wind termination shock 
to agree with the afterglow shock radius inferred from the time of the first 
bump. When p = 2 the density profile rises by almost an order of magnitude 
and then decreases, more gradually, back to its initial value. The initial rise 
agrees with the one suggested by Lazzati et al. (2002), however, Lazzati et al. 
suggest a consequent decrease in the density to a factor of 5 below the initial 
density value followed by a second and smaller density bump, where according 
to our results such a large dip in the density is not required. The difference 
between the profiles arises mainly due to the different approximation used for 
the angular smoothing effect (see the Appendix). 

So far, we have assumed a spherically symmetric external density profile, n = 
n(r). This may occur due to a variable stellar wind, but is not expected for 
an ISM. As we obtain that an underlying constant density profile provides a 
better fit for GRB 0210004, it is more natural to expect density fluctuations 
in the form of clumps, rather than being spherically symmetric, in this case. 
This interpretation requires p ~ 2 for which the density at large radii decreases 
back to its value at i?o- As the density profile for p = 2 is not smooth, several 
density clumps are needed. The first clump should be at Ri ~ 1.5R and 
with an over-density of factor ~ 8. In order to have a similar effect as a 
spherical density bump, the clump must replace all the emitting material, 
i.e. its size (radius), Z c i, must be large enough so that its mass is larger than 
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Fig. 1. The external density profiles, n(R), that reproduce the R-band light curve, 
for different values of electrons' index, p: p = 2.4 (dashed-doted), p = 2.2 (dashed 
line) and p = 2 (solid line). The thick lines are the exact numeric solution of Eqs. 
8-10. The thin lines are the analytic solution of Eq. 12. The thin dotted line depicts 
the expected amplitude of the density fluctuations, An(R) oc R~ 5 ^ 2 , for a random 
distribution of clumps inside a uniform density background, no. The normalization is 
derived assuming that the first bump is due to a single clump. The inset on the right 
depicts the R-band observed data points and the fitted light curve. The observed 
R-band data points are taken from: Fox 2002, Uemura et. al. 2002, Oksanen k Aho 
2002, Rhoads et. al. 2002, Winn et. al. 2002, Zharikov et. al. 2002, Halpern et.al. 
2002a k 2002b, Balman et. al. 2002, Cool k Schaefer 2002, Holland et. al. 2002a k 
2002b, Bersier et. al. 2002, Sahu et. al. 2002, Oksanen et. al. 2003, Matsumoto et.al. 
2002, Stanek et. al. 2002, Mirabel et. al. 2002, Masseti et. al. 2002, Barsukova et. al. 
2002 , Malesani et. al. 2002, Mirabel et. al. 2002. 

the swept up mass at that radius within an angle of 7 _1 around the line of 
sight: l c \ > Z min = (n /n c i47 2 ) 1/3 _Ri « 0.03i?i « 10 16 cm. An upper limit on 
the size of the clump can be put from the fact that the bump in the light 
curve decays on a time scale At ~ t. Since R oc t 1 / 4 for an ISM, this implies 
l c \ < ^max = (2 1/4 - l)Ri « 0.19.Ri (Lazzati et al. (2002) obtain a similar 
clump size using different considerations). 

Assuming a homogeneous distribution of clumps with the same physical size 
and over-density, the mean distance between neighboring clumps is d c \ ~ 
(ttI^Ri) 1 / 3 « 4 — 5 x 10 16 cm, where the numerical estimate assumes l c \ = 
in which case the clumps hold roughly 5% of the volume and 30% of the mass 
(these are lower limits as l c \ > l m i n would imply larger filling factors). There- 
fore, soon after the collision with the first clump we expect overlap between 
pulses from different clumps, where the number of clumps that intersect a 
given shell with a radius R and angular size I/7 is N c \ ~ R 2 l c \/^ 2 d^ x oc R 5 . 
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Since, on average, the clumps hold a constant fraction of the shell's mass, a 
single clump constitute a fraction oc R~ 5 of the matter at this radius. The to- 
tal fluctuation in the density would be, therefore, oc R~ 5 ^ 2 . This is in a rough 
agreement with the fluctuations in the density profile we have obtained for 
p=2 (see Fig 1). 



3.2 A Variable Energy Profile 

We solve Eqs. 14 and 15 numerically, for a constant ISM density profile, as- 
suming that the energy is constant, E , up to the first observation at to, and 
letting E vary from this point onwards. Figure 2 depicts the energy profile 
obtained for different values of p as a function of 9 = I/7 (the angular size 
of the observed area). An electron power law index of p = 2.6, requires an 
almost monotonous increase of E in the observed region. Such a profile may 
arise due to refreshed shocks. However, the continues increase in E requires a 
continuous arrival of new shells, a scenario which we consider as unlikely. The 
energy profiles obtained for p = 2.2 and p = 2.4 could reflect irregular patches 
with an initial angular size of 9q ~ 0.02 rad and an average energy of sev- 
eral times Eq. The energy fluctuations decrease with time, as expected from a 
patchy shell (see Figure 2) . The profile obtained for p = 2 shows an initial rise 
followed by a gradual (and bumpy) decrease back to the initial value. Such a 
profile can correspond to a line of sight is ~ 0.04 rad away from a hot spot 
(the average energy over a large area is E ). This hot spot may be a hot patch 
in an irregular jet. Alternatively, as suggested by Lazzati et al. (2002), this 
hot spot may be the core of a jet (on the jet axis) in an axisymmetric angle 
dependent regular jet 7 . According to this interpretation the angular size of 
the jet's core is 9 C ~ 0.02 rad, the isotropic equivalent energy outside the core 
is roughly constant and its value is ~ 3 times less than the core's energy. 



4 Discussion 

We have presented general expressions for the afterglow light curve when the 
energy in the blast wave varies with time and for a variable external density 
profile. This formalism follows and generalizes the work of SPN98, and relates 
the variability in the energy and density to the variability in the light curve. 
Despite the variability in the light curve, the shape of the broad band spectrum 



7 Though the wiggles in E{9) require some additional small amplitude variability 
on small angular scales on top of an underlying smooth axisymmetric jet profile on 
large angular scales. 
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Fig. 2. The isotropic equivalent energy, E, within an angle 9 = I/7 around the line 
of sight as a function of 6, for different values of p: p = 2 (solid line), p = 2.2 (dashed 
line), p = 2.4 (dashed-dotted line) and p = 2.6 (bold dotted line). The curves are 
normalized by the value at the first observation, Eq. The thin dotted lines outline the 
expected fluctuations in E for a patchy shell model with fluctuations on an angular 
scale fl = 0.02 (using E = 2AE ± 1.9E o (0 & /9)). The inset on the left depicts the 
observed data points as a function of 6 and the fitted curve. 

remains the same, with some variability in the values of the break frequencies 
and flux normalization. 

We have focused on the slow cooling spectrum at frequencies v > z/ m , and 
derived detailed equations for these they seem the most relevant for 

the majority of observed optical light curves. Similar equations can be easily 
derived for other spectral regimes using Eqs. 4-6. We find that for v m < v < u c , 
variability in the light curve can be induced both by variability in the energy or 
by variability in the external density (or both) . A similar behavior is expected 
for v < min(z/ m , u c ), for both slow and fast cooling. For v > max(z/ m , v c ) we 
find that a variable density hardly induces any fluctuations in the light curve, 
while a variable energy can induce significant fluctuations. We expect a similar 
behavior for u c < u, in the fast cooling regime. 

We applied our formalism to GRB 021004, which displayed significant devia- 
tions from a simple power law decay in its optical (R-band) light curve. We find 
that several different models may provide a reasonable fit to the observed light 
curve. These include models where the variability is induced either by density 
fluctuations or by energy fluctuations, where the latter may be caused either 
by refreshed shocks or by a patchy angular structure of the GRB outflow. 
These models vary significantly with the value of p. Chandra's observations 
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t days 

Fig. 3. The expected light curve F v {t)/F v fl for v > v c where F u> q is the observed flux 
at to = 537 sec. The expected light curves are calculated using (i) the energy profile 
that reproduces the R-band light curve for p = 2 (solid line) and (ii) the density 
profile that reproduce the R-band light curve for p = 2 (dashed line). The short thick 
line represents Chnadra's X-ray measurement, normalized for the expected flux at 
20.5 hr in the varying energy light curve (a power law decay with index —1 from 
20.5 hr till 44 hr after the trigger). 

constrain the electron's index to be p = 2.2 ± 0.2, but even under this con- 
strain many different models can produce the observed light curve. A tighter 
constrain would limit the models considerably. The following models provide 
a viable fit to the light curve: I) A variable density: a) For p = 2.2 there is an 
order of magnitude rise in the density followed by a roughly constant density; 
b) For p = 2 we find a similar rise, but then the density gradually decreases 
back to its initial value; II) A variable energy: a) For p = 2.6 refreshed shocks 
are required in order to explain the energy profile; b) For p ~ 2.2 — 2.4 a 
patchy shell model provides a good fit; c) For p = 2 a hot spot (possibly the 
core of an axisymmetric jet) should reside near our line of sight. 

As any given single band light curve (which does not show a strong variability 
on time scales At <C t) can be reproduced by either density or energy varia- 
tions, it is important to find ways to distinguish between these two models and 
their variants. An independently determined value of p, say from the spectrum, 
would have made this task easier (but still not completely determined). Simul- 
taneous light curves both above and below the cooling frequency, is c , provide 
the best way to differentiate between a variable energy and a variable density: 
for the latter strong variability is possible only below v c . Figure 3 depicts the 
light curves that are predicted above v c , using the energy or density profiles 
deduced from the R-band light curve, that is assumed to be below v c . Chan- 
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dra obtained an X-ray light curve between 1 and 2 days (Sako et al. 2002, the 
thick lines in Figure 3). Unfortunately, by this time the fluctuations expected 
in the X-ray light curve according to the two models are rather similar and 
it is hard to distinguish between them. Still, it would be interesting to search 
for a correlation between the R-band light curve and the X-ray light curve at 
this time. An earlier X-ray observation could have enabled a clear distinction 
between the two models. 

A variable energy model could arise either from refreshed shocks or from angu- 
lar inhomogeneity in the jet. In the refreshed socks scenario, we expect during 
the collision between the two shells an increase in the spectral slope (3 (defined 
by F v oc v 13 ) and a strong signal in the radio, (Kumar & Piran 2000a). This 
emission should last over At ~ t. A refreshed shocks can only add energy to 
the blast wave the total energy in this picture can only increase with time. 
In the patchy shell model we expect random fluctuations whose amplitude 
decays with time as I/7 (see Fig. 2). 

Although the current observations do not enable us to determine which one 
of the scenarios described above is the correct one (if any), we feel that the 
patchy shell model with p = 2.2 (which agrees with the p = 2.2 ± 0.2 value 
suggested by Chandra's observations) is the most likely scenario. According to 
this interpretation the line of sight of GRB 021004 falls in a "cold spot" where 
the energy is 2.5 times below the average. This agrees with the observation of 
rather low 7-ray flux from this burst. The total 7-ray energy, — 3 • 10 50 ergs 
is within the standard deviation of the energy distribution presented by Frail 
et. al (2001), but it is 1.5 times smaller than the average value. On the other 
hand an extrapolation of Chandra's measured X-ray flux (Sako et al. 2002) 
to 11 hr after the burst yields F x pa 9 • 10" 13 ergs cm -2 sec -1 . This value 
is 1.5 times larger than the narrowly clustered value of F x in other bursts: 
6 • 10" 13 ergs cm" 2 sec" 1 (Piran et. al. 2001). The X-ray flux reflects E k , the 
kinetic energy of the relativistic ejecta (averaged over an angular scale I/7 
corresponding to 7 ~ 10). Hence, in this burst E k /E^ is larger by a factor of 
2.25 than the average value. This factor is similar to the energy fluctuations 
we find in the patchy shell model for p = 2.2 (see Fig. 2). While in most GRBs 
that show a larger value of E k /E 1 we, most likely, observe a 7-rays hot spot 
(Piran 2001). According to this interpretation GRB 021004 is the first burst 
in which a clear 7-ray cold spot has been seen. 
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This work was partially supported by the Horwitz foundation (EN) and by 
the Institute for Advanced Study, funds for natural sciences (JG). 

Appendix 

Because of the curvature of the afterglow shock, that is spherical rather than 
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planar, photons that are emitted from the shock front at the same time in 
the source rest frame (i.e. at the same radius), but at different angles from 
the line of sight, reach the observer at different times. This causes two main 
effects: first, the bulk of the energy that is emitted at a given time at the 
source is delayed compared to a photon emitted on the line of sight at that 
time, and second, at any given time the observer receives photons that were 
emitted at different radii. In our analysis we take the first effect (angular time 
delay) into account (see Eq. 3), but the second effect (angular smoothing) is 
neglected (see Nakar & Piran 2003 for a solution of the spherical symmetric 
afterglow light curve that takes a full account of the angular effects). For 
spherical shells, the angular smoothing produces an observed light curve which 
is a smoothed version of the line of sight emission. The relative importance of 
angular smoothing is determined by the ratio t\ os /t ang , where U os = dr/'-f 2 
and t ang ~ R/4cy 2 . When the external density decays as a power law, n oc r~ k , 
the line of sight time is: t\ os = R/4(4 — k)cry 2 . Most of the contribution to the 
observed flux at a time t Q b s , comes from emission at radii 4 < R/ ct Q ^ 2 {R) < 
4(4 — k), which correspond to t\ os < t Q b s < tang- Hence, this effect is important 
when the light curve from the line of sight varies significantly (compared to 
the smooth power law decay) on time scales shorter than t ang (i.e. Ati os /ti os < 
^angAios ~ 4), which corresponds to density variations on AR/R < 0.4. In 
such the observed light curve is significantly less variable than the line 

of sight light curve. 

We calculate the density profile assuming that the observed (smoothed) light 
curve is similar to that from the line of sight. Thus, the real density profile 
has to be more variable than the one we obtain. The difference between the 
two is smaller when the density profile increases with radius, and the emission 
along the line of sight increases with time (compared to the power law decay) . 
In this case, the observed flux is dominated by emission from large radii, near 
the line of sight (with a relatively small contribution from large angles, for 
which the emission took place at smaller radii, where the external density was 
relatively low) and the angular smoothing effect is less important. However, 
when the density drops, the angular effect becomes important. Panaitescu & 
Kumar (2000) have shown that even a sharp drop in the density produces only 
a gradual temporal decay in the observed light curve, and that the angular 
smoothing dictates a maximal power law index of the temporal decay at late 
times. 

The observed R-band light curve of GRB 021004 shows variations on time 
scales of At ft ~ 1 < £ a ngAios ~ 4, therefore the angular smoothing effect 
is not negligible. This effect can be seen by comparing our density profile to 
the density profile obtained by Lazzati et al. 2002 (their Fig. 1), which take 
the angular smoothing effect into account. The main difference between the 
profiles is in the sharp density drop after the first density bump (the density 
in the profile of Lazzati et al. drops to one order of magnitude below the initial 
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density, while our profile drops back to the initial density). This sharp drop 
is required, when the angular smoothing is considered, in order to obtain the 
steep temporal decay in the light curve after the first bump (t ~ 10 4 sec). 

This last result is obtained under the assumption of a full spherical symme- 
try. However, it is more likely that the overdense regions are concentrated in 
clumps and not in spherical shells (see section 3.1 and Lazzati et al. 2002). 
The radial size of the first clump is ~ 10 16 cm ~ Ri/li, where Ri and 7$ are 
the radius and Lorentz factor of the blast-wave when it first interacts with 
the clump. If we assume that the clump is spherical then its angular size is 
~ 1/7; and at the beginning of the first bump (t ~ 10 3 sec) it "fills" most 
of the observed region (the region within an angle of I/7 around the line of 
sight). As the dominant emission is from the hot clump, the angular time 
remains constant, R/Ajf, while t\ os grows. Therefore, the angular smoothing 
effect becomes less and less important, and the approximation which neglects 
this effect holds better than for a spherically symmetric external density pro- 
file. Therefore, in this scenario, our method yields a good approximation for 
the actual density profile at radius R, averaged over an angle of ~ I/7 around 
the line of sight. 
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